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Abstract 

Irreversible opinion spreading phenomena are studied on small- 
world and scale-free networks by means of the magnetic Eden model, 
a nonequilibrium kinetic model for the growth of binary mixtures in 
contact with a thermal bath. In this model, the opinion of an indi- 
vidual is affected by those of their acquaintances, but opinion changes 
(analogous to spin flips in an Ising-like model) are not allowed. We 
focus on the influence of advertising, which is represented by external 
magnetic fields. The interplay and competition between temperature 
and fields lead to order-disorder transitions, which are found to also 
depend on the link density and the topology of the complex network 
substrate. The effects of advertising campaigns with variable dura- 
tion, as well as the best cost-effective strategies to achieve consensus 
within different scenarios, are also discussed. 

Keywords: Complex networks; sociophysics; kinetic growth models; 
phase transitions; nonequilibrium processes. 



1 Introduction 

Many interdisciplinary fields of science are nowadays benefitting from valu- 
able tools, insight and procedures borrowed from long-standing, traditional 
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applications of statistical physics [H [21 El 111 E] ■ In particular, a rich variety of 
social phenomena, such as social influence and self-organization, cooperation, 
opinion formation and spreading, have recently been studied using models 
and methods of statistical physics (see e.g. [HllTlElElliniiniinilinilTlIISl 
fT6l [T71 fT8l fT9| EOl EI])- In some of these investigations, social processes were 
studied by means of Ising-type spin models [221l231[21[2Sl[2Sll2Zl[2Hll2n],in 
which spin states denote different opinions or preferences, while the coupling 
constant describes the convincing power between interacting individuals in 
competition with the "free will" given by the thermal noise. 

Four decades ago, seminal observations by Milgram [30] showed that the 
mean social distance between a pair of randomly selected individuals was 
astonishingly short, typically of a few degrees. Since then, the phenomenon 
of average shortest path-lengths that scale logarithmically with the system 
size, known as small-world effect, has been found ubiquitously in many dif- 
ferent real networks, which also generally show highly interconnected local 
neighborhoods. 

The well-studied classical random graphs, which are networks built by 
linking nodes at random, display the small-world effect but have much lower 
connectivities than those usually observed in real networks. To address this 
problem, the so-called small-world networks were proposed few years ago [311 
[32] as a realization of complex networks having short mean path-lengths (and 
hence showing the small-world effect) as well as large connectivities. Starting 
from a regular lattice, a small-world network is built by randomly adding or 
rewiring a fraction p of the initial number of links. Even a small fraction of 
added or rewired links provides the shortcuts needed to produce the small- 
world effect, thus displaying a global behavior close to that of a random 
graph, while preserving locally the ordered, highly connected structure of a 
regular lattice. 

Small-world networks, showing the appropriate topological features ob- 
served in real social networks, can thus be meaningfully used as spin model 
substrates in order to investigate different social phenomena. In this context, 
irreversible opinion spreading phenomena on small- world [331 [Mj and other 
complex networks [33] have been studied by means of the so-called magnetic 
Eden model (MEM), a nonequilibrium kinetic growth model in which the 
deposited particles have an intrinsic spin and grow in contact with a thermal 
bath [371 [3H] (for a review, see Ref. [39]). The MEM defined on small- 
world networks present continuous order- disorder phase transitions taking 
place at finite critical temperatures for any value of the rewiring probability 
p > [32], a phenomenon analogous to observations reported previously in 
the investigation of equilibrium spin systems [231 [21] • 

According to the growth rules of the MEM, which are given in the next 
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Section, the opinion or decision of an individual are affected by those of 
their acquaintances, but opinion changes (analogous to spin flips in the Ising 
model) do not occur. It should be noticed that this behavior contrasts with 
other models of opinion formation and evolution, which generally consider the 
possibility of changing opinion. As such, the MEM is applicable to scenarios 
were opinion thermalization processes are unlikely to occur. 

Within this context, the aim of this work is to focus on the influence of 
advertising, which is represented by an external, homogeneous magnetic field. 
The interplay and competition between temperature and fields lead to order- 
disorder transitions, which are investigated by means of order-parameter 
probability distributions and ensemble averages. The corresponding tem- 
perature vs. magnetic field order- disorder phase diagram is obtained, and 
the role played by shortcuts is discussed. Moreover, the effects of advertising 
campaigns with variable duration, as well as the best cost-effective strategies 
to achieve consensus within different scenarios, are also investigated. 

Furthermore, in order to explore effects arising from the nature of the 
complex network substrate, we examine the MEM growing on scale-free net- 
works. The degree distribution of these networks, P{k), which measures the 
probability that any randomly chosen node has k links, follows a power-law, 
i.e. P{k) ~ k~'^. Unlike the Poisson degree distribution of classical random 
graphs and small-world networks, power-law distributions lack any charac- 
teristic scales. The structure of scale-free networks provides an interesting 
setting in which social roles differ significantly among individuals. Indeed, 
the most connected nodes in the scale-free network (usually designated as 
hubs) represent highly influential social leaders, whose opinions can poten- 
tially affect the decisions of many other individuals in the society. 

This paper is organized as follows: in Section 2, details on the model 
definition and the simulation method are given; Section 3 is devoted to the 
presentation and discussion of the results, and Section 4 contains the conclu- 
sions. 

2 Model and simulation method 

In this work, we consider two paradigmatic types of complex network model: 
(a) the one-dimensional, nearest-neighbor, adding-type small-world network 
(SWN) IIQIIIT], and (b) the Barabasi- Albert (BA) scale-free network (SF) 

1121 1131113. 

In case (a), starting with a ring of sites and A^ bonds, a SWN realiza- 
tion is built by adding new links connecting pairs of randomly chosen sites. 
For each bond in the original lattice, a shortcut is added with probability 
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p. During this process, multiple connections between any pair of sites are 
avoided, as well as connections of a site to itself. Since the original lattice 
bonds are not rewired, the resulting network remains always connected in 
a single component. On average, pN shortcuts are added and the mean 
coordination number is (z) = 2(1 + p). 

In case (b), SF networks are built by following the preferential attachment 
growth mechanism: starting with a small fully-connected graph of size m, 
a new node with m edges is added at every time step and is connected 
to m different nodes already present in the system, where the probability 
for an already existing node to acquire a new link is proportional to its 
degree. The BA type of scale-free network is characterized by power-law 
degree distributions, i.e. P{k) ~ k~^, where the degree exponent is 7 = 3. 

Once the network is created, a randomly chosen up or down spin is de- 
posited on a random site. Then, growth takes place by adding, one by one, 
further spins to the immediate neighborhood (the perimeter) of the growing 
cluster, while taking into account the corresponding interaction energies. By 
analogy to the Ising model, the energy of a configuration of spins is given 



where Si = ±1 indicates the orientation of the spin for each occupied site 
(labeled by the subindex i), J > is the ferromagnetic coupling constant 
between nearest-neighbor (NN) spins, > is the external field, and {ij) 
indicates that the summation is taken over all pairs of occupied NN sites. As 
with other spin systems defined on complex networks, the magnetic interac- 
tion between any pair of spins is only present when a network link (either a 
lattice bond or a shortcut) connects their sites. 

Setting the Boltzmann constant equal to unity {ks = 1), the probability 
for a new spin to be added to the (already grown) cluster is defined as 
proportional to the Boltzmann factor exp{—AE/T), where AE is the total 
energy change involved and T is the absolute temperature of the thermal 
bath. Energy, magnetic field and temperature are measured in units of the 
NN coupling constant, J, throughout. At each step, all perimeter sites have 
to be considered and the probabilities of adding a new (either up or down) 
spin to each site must be evaluated. Using the Monte Carlo simulation 
method, all growth probabilities are first computed and normalized, and then 
the growing site and the orientation of the new spin are both determined by 
means of a pseudo-random number. Although the configuration energy of 
a MEM cluster, given by Eq.([l]), resembles the Ising Hamiltonian, it should 
be noticed that the MEM is a nonequilibrium model in which new spins are 
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continuously added, while older spins remain frozen and are not allowed to 
flip. The growth process naturally stops after the deposition of N particles, 
when the network becomes completely filled. 

Since finite- size effects for the MEM growing on small- world and BA 
scale-free networks are mild (see Refs. [33l |35] for detailed discussions), in 
this work we adopt a fixed network size, N = 10^, throughout. For any given 
set of defining parameters (i.e. the shortcut-adding probability p, the BA 
parameter m, the magnetic field H and the temperature T), ensemble aver- 
ages were calculated over 10^ different (randomly generated) networks, and 
considering typically 50 different (randomly chosen) seeds for each network 
configuration. 

3 Results and discussion 

3.1 Advertising and opinion spreading in small-world 
networks 

The natural order parameter of a magnetic system is the total magnetization 
per site, i.e. 



which, in the context of this work, is to be measured on the completely 
filled network. The normalized probability distribution of the magnetization, 
Pm, is shown in Figure 1 for different temperatures, shortcut densities and 
magnetic fields: (a) p = and H = 0.05; (b) p = and H = 0.2; (c) p = 0.5 
and H = 0.05; (d) p = 0.5 and H = 0.2. 

For the regular lattice (Figures l(a)-(b)), the distributions undergo abrupt 
changes when passing from low to high temperatures. In the case of small 
fields (Figure 1(a)), the low-temperature distribution is dominated by sharp 
maxima taking place at M = ±1. Increasing temperature, these peaks re- 
tract while a local maximum develops near M ^ 0.2 forming a Gaussian-like 
distribution. Notice also that, due to the small but positive magnetic field, 
the distributions are skewed and biased towards positive values of the mag- 
netization. For large fields (Figure 1(b)), the distributions are highly asym- 
metric and show the same kind of abrupt transition from a low-temperature 
sharp maximum at M = 1 to high-temperature Gaussian-like shapes. 

"'^ Notice that, although this work is mainly motivated by social phenomena, a magnetic 
language is adopted throughout. Indeed, physical concepts such as temperature and mag- 
netization, spin growth and clustering, ferromagnetic-paramagnetic phase transitions, etc, 
can be meaningfully re-interpreted in sociological/sociophysical contexts. 




(2) 
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Figure 1: Normalized probability distributions of the magnetization for dif- 
ferent temperatures, as indicated. The panels correspond to different short- 
cut densities and magnetic fields: (a) p = and H = 0.05; (h) p = and 
H = 0.2; (c) p = 0.5 and H = 0.05; (d) p = 0.5 and H = 0.2. Sharp maxima 
at M = ±1 have been truncated in panels (a)-(b). 



In contrast, for the small- world network under small fields (Figure 1(c)) 
one observes the onset of two maxima located at M = — M„eg and M = Mpos 
(such that < M^eg < Mpos < 1)5 which become sharper and approach 
M = ±1 as T is gradually decreased. The smooth shift of the distribution 
maxima across T c:^Tc, from M at high T to the low-temperature nonzero 
spontaneous magnetizations M = —Mn^g, Mpos, is indeed the signature of 
true thermally-driven continuous phase transitions [3H1I1S]- Hence, the pres- 
ence of shortcuts in the small-world network triggers a critical behavior in 
the irreversible growth of MEM clusters [33J , analogously to phenomena ob- 
served in related equilibrium systems such as the Ising model [231 [2^ . When 
strong fields are applied (Figure 1(d)), the negative maxima at M = —M^eg 
vanish, thus leading to highly asymmetric distributions peaked at M = Mpos 
that tend to M„„, ~ 1 as T decreases. 
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Figure 2: Ensemble-averaged absolute magnetization as a function of the 
magnetic field for different shortcut densities and temperatures, as indicated. 

Since, for if ^ 1, the low-temperature probability distributions of the 
magnetization are nearly symmetric and sharply peaked around M ^ ±1, 
the ensemble-averaged total magnetization vanishes. This well-known short- 
coming due to finite-size effects is avoided by considering the absolute value 
of the magnetization, |M|, as the actual order parameter |l5l 1^ . 

Figure 2 shows the ensemble-averaged absolute magnetization as a func- 
tion of the applied magnetic field for different shortcut densities and temper- 
atures, as indicated. For any value of the temperature, the order parameter 
plot for the small-world network (solid lines) lies above the corresponding 
one for the lattice (dotted/dashed lines), since shortcuts favor long-range or- 
dering connections between distant clusters across the network. Moreover, 
as expected, most of the plots are monotonically increasing due to the or- 
dering effect of the external field. Interestingly, however, a slight decrease in 
the absolute magnetization is observed at intermediate fields, H = 0.1 — 0.3, 
for p = 0.5 and T = 0.3. This behavior is indeed due to the competition 
between the strong clustering tendency of neighboring spins at low tempera- 
tures, which favors the formation of extended parallel-oriented up and down 
domains, and the external magnetic field. 

Figure 3 shows p— dependent, order- disorder transition curves on a H 
vs T phase diagram, where transition states are defined by the condition 
(|M|) = 0.5. Due to the long-range ordering effect of shortcut connections, 
transitions for small- world networks with higher shortcut densities are shifted 
towards higher temperatures. Moreover, due to the global reach of the ex- 
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Figure 3: Phase diagram in the H vs T plane for different values of the 
shortcut fraction p, as indicated. 

ternal field, the ordered phase is attainable even for very high temperatures. 
According to the sociophysical interpretation of this model, we conclude that 

consensus can be reached irrespective of the degree of intrinsic cohesion and 
intcrconnectedness (as given by the parameters T and p, respectively) as long 
as a sufficiently intense, global mass media advertising campaign is carried 
out. 

Figure 4 shows the transition temperature as a function of the shortcut 
fraction p for different values of the external field, as indicated. The mono- 
tonic shift towards higher temperatures for small- world networks with higher 
connectivities is observed to be roughly independent of the external field. In 
order to exploit this approximate scaling property, we define an effective, 
shortcut- induced temperature shift ATe//(p) — —ap and re-plot the phase 
diagram in terms of H vs T^ff = T + AT^jf, which is shown in the inset to 
Figure 4. Using a = 0.45, the small-world network data tend to collapse into 
a curve that can be well fitted by the relation H = 0.27 x T^ff (solid line). 
As expected, lattice data show a distinct behavior: the p — transition is 
well fitted by if = 0.32 x T^ff (dashed line), which clearly departs from the 
small-world network data collapse. 

In order to gain further insight into the dynamical features of the opinion 
spreading process, let us consider the effects of advertising campaigns with 
variable duration. For the sake of simplicity, we assume a constant deposition 
(or spreading) rate and measure time ast = n/N, i.e. the fraction of already 
deposited particles, n, relative to N, the system's final size. At every stage 
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Figure 4: Order-disorder transition temperature as a function of the shortcut 
fraction p for different values of the external field, as indicated. The inset 
shows the phase diagram H vs Tf.ff, where the effective temperature Tg// — 
T — 0.45 X p leads to the collapse of the small-world network data. Separate 
fits to the lattice (dashed line) and small-world network (solid line) data are 
also shown. 

during the growth process, we determine whether a given configuration has 
a majority of up spins {uup > n/2) and, after ensemble- averaging over many 
configurations, we measure the fraction of states with a majority of up spins, 
fup- Since the seeds are chosen at random, the initial condition is /„p = 0.5 
and the success of the campaign can be measured by the extent to which 
fup approaches unity. We assume that the field is applied from t — Q to 
a specified time t — tn] afterwards, the field is switched off. Hence, the 
parameter < < 1 represents the duration of the advertising campaign. 

Figure 5 shows time profiles of f^p for small-world networks with p = 0.5, 
different temperatures and magnetic fields: (a) T = 0.5 and H = 0.1; (b) 
T = 0.5 and H = 0.5; (c) T = 2 and H = 0.1; (d) T = 2 and H = 0.5. 
Different field switch-off times correspond to different plots on each panel, 
as indicated. Figure 5(a) shows that small fields are not capable of affecting 
the strong clustering tendency typical of low temperatures. Indeed, the mag- 
netization probability distributions are sharply peaked at Mneg ~ Mpos ~ 1 
(similarly to the low-temperature plot in Figure 1(c)) and are little infiuenced 
by the small external field even in the case oi tn — For low temperatures 
and strong fields, instead, fup grows steadily with tn (Figure 5(b)), since 
strong fields deplete the probability distribution for negative values of the 
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Figure 5: Time profiles for tlie fraction of states with a majority of up spins, 
fup, for small- world networks with p — 0.5 and different field switch-off times, 
tn, as indicated. The panels correspond to different temperatures and mag- 
netic fields: (a) T = 0.5 and H = 0.1; (b) T = 0.5 and H = 0.5; (c) T = 2 
and H = 0.1; (d) T - 2 and H = 0.5. 



magnetization (recall Figure 1(d)). Moreover, f^p remains constant after the 
field is switched off due to the lack of sizable thermal fluctuations. Notice 
that highly successful final results (i.e. fup ~ 1 for i — > 1) can be achieved 
with intense but very short campaigns {tn < 0.1). Figure 5(c) shows that, 
for high temperatures and small fields, fup strongly depends on the cam- 
paign duration: /„p increases monotonically with tn and is close to unity for 
t^f = 1. Comparing to Figure 5(a), one observes that higher temperatures 
increase the overall response sensitivity of the system under the infiuence of 
small fields. The case of high temperatures and large fields is shown in Fig- 
ure 5(d). The strong fields naturally drive the system towards large values of 
fup] however, whenever the fields are switched off, the trend towards thermal 
disorder sets back on and fup decreases again. 

According to these results, the best cost-effective strategy to achieve a de- 
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Figure 6: Dependence of the ensemble-averaged final states [t = 1) on the 
field switch-off time, tn, for different temperatures, external fields and ob- 
servables: (a) the total magnetization, (b) the absolute magnetization, and 
(c) the fraction of final states with a majority of up spins. These results are 
for small- world networks with p = 0.5. 



sired majority consensus {fup — 1 for t — > 1) depends crucially on the degree 
of intrinsic cohesion tendency of the society, here modeled by the tempera- 
ture. In the case of strongly cohesive social groups (low T), intense, short 
campaigns yield larger consensus levels than mild, long ones. Interestingly, 
however, the opposite is true in the scenario of weakly cohesive societies (high 
Tg 

^Notice that the low- and high-temperature regions are relative to the corresponding 
zero-field critical temperature. For SWNs with N ~ 10"^ and p = 0.5, the zero-field 
transition is located at Tc « 0.6 I33j. Also recall that temperatures are measured in units 
of the nearest-neighbor coupling constant throughout (see Sect. 2 for details). 
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Figure 7: Ratio fup{T = 0.5)//„p(T = 2) as a function of the field switch-off 
time tn for small- world networks with p = 0.5 under different external fields, 
as indicated. 



Figure 6 shows the dependence of the ensemble-averaged final states (i.e. 
t = 1) on the campaign duration (or, equivalently, the field switch-off time, 
tn) for different temperatures, external fields and observables: (a) the total 
magnetization, (b) the absolute magnetization, and (c) the fraction of final 
states with a majority of up spins. The results shown correspond to small- 
world networks with shortcut fraction p = 0.5. 

For low temperatures and strong fields (dash-dotted line plots), the sys- 
tem becomes rapidly ordered in the up direction and stays in the up-ordered 
state even after the field is withdrawn. Hence, (M) and fup follow a similar 
pattern, namely a monotonic, steadily increasing behavior as tn is increased. 
The absolute magnetization is close to unity but displays a dip, which is 
due to the passing of the down-ordered probability distribution, through the 
zero-magnetization region, towards the up-ordered side. This behavior is 
similar to the dip previously observed in the thick solid line plot of Figure 
2. For low temperatures and small fields (dotted lines), the positive fields 
are not capable of significantly affecting the strongly clustered down-ordered 
domains, thus leading to small but moderately increasing values for (M) 
and fup- In contrast, (|M|) is large but slightly decreasing, as explained by 
mild shifting effects of the magnetization probability distribution towards the 
positive side. When high temperatures are considered (long-dashed and solid 
lines), the magnetization probability distributions are nearly Gaussian with 
positive mean values. Hence, (M) and (|M|) are expected to be similar and 
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Figure 8: Phase diagram in the H vs T plane for different values of the 
Barabasi- Albert parameter m, as indicated. 

generally small. However, since in the — > limit these distributions are 

Gaussians centered at M = 0, even slight shifts towards M > when small 
fields are applied lead ultimately to fairly large values of This explains 
the fact that, for high temperatures, mild, long campaigns yield higher f^p 
values than short, intense ones. 

Figure 7 shows the ratio fup{T — 0.5)//„p(T = 2) as a function of tn for 
different fields. For H = 0.1, the ratio remains fairly constant about unity 
up to tj:/ ~ 3 X 10~^, while it decreases monotonically for higher values of 
tn and reaches a minimum aX tn = I (horizontal arrow), where the high- 
temperature campaign is more successful by more than 40% relative to the 
low-temperature one. Hence, mild, long campaigns are optimal for (high- 
temperature) weakly cohesive societies. For H — 0.5, the ratio peaks at tn — 
3 X 10~^ (vertical arrow) and then decreases towards unity for larger values of 
tn, thus showing that for (low-tcmpcrature) highly cohesive societies, short, 
intense campaigns are the best strategy. 

3.2 Advertising and opinion spreading in scale-free net- 
works 

Unlike the case of small-world networks, the highly inhomogeneous degree 
distribution of scale- free networks represents a scenario in which social roles 
differ significantly among individuals. In particular, hubs in the scale-free 
network represent highly infiuential social leaders, whose opinions can po- 
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Figure 9: Order-disorder transition temperature as a function of the 
Barabasi-Albert parameter m for different values of the external field, as 
indicated. The inset shows the phase diagram H vs Tgff, where the effective 
temperature T^jf = T — 0.50 x m leads to an approximate data collapse for 
scale-free networks with m > 1. 



tentially affect the decisions of many other individuals in the society. 

The paradigmatic type of scale-free network is based on the preferential 
attachment growth mechanism: a new node with m edges is added at every 
time step and is connected to m different nodes already present in the system, 
where the probability for an already existing node to acquire a new link is 
proportional to its degree. This simple rule leads to an ensemble of Barabasi- 
Albert (BA) scale-free networks, which are characterized by power-law degree 
distributions with exponent 7 = 3. 

Figure 8 shows order-disorder transition plots on a i7 vs T phase diagram 
for the MEM growing on SF networks with different values of the BA param- 
eter m. Since the overall average degree grows linearly with m, the transition 
curves for SF networks with higher hnk densities are shifted towards higher 
temperatures. Moreover, due to the global reach of the external field, the 
ordered phase is attainable even for very high temperatures. Comparing to 
the SWN results (see Figure 3), we observe that, irrespective of the details 
of the underlying complex network substrate, consensus and order can be 
reached for any level of intrinsic cohesion and interconnectedness. 

Figure 9 shows the transition temperature as a function of the BA pa- 
rameter m for different values of the external field, as indicated. Except for 
the m — 1 case, the monotonic shift towards higher temperatures appears 
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Figure 10: Dependence of the ensemble-averaged final states [t = 1) on 
the field switch-off time, t^, for different temperatures, external fields and 
observables: (a) the total magnetization, (b) the absolute magnetization, and 
(c) the fraction of final states with a majority of up spins. These results are 
for scale-free networks with BA parameter m = 3. 



to be roughly independent of the external field, which suggests to define an 
effective, m-induced temperature shift ATe//(m) = —am. Using a = 0.5, we 
obtain an approximate data collapse for scale- free networks with m > 1, as 
shown in the inset to Figure 8. Comparing to the small-world network data 
shown in Figure 4, we confirm that, from a qualitative point of view, the 
critical behavior induced by random connections is generally characteristic 
of complex network substrates. 

Figure 10 shows the dependence of the ensemble-averaged final states 
as a function of the field switch-off time, t^, for different temperature^. 



^Since, for BA-SF networks with N = 10^ and m — 3, the zero-field transition is located 
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Figure 11: Ratio fupiT = 1)/ fup{T = 4) as a function of the field switch-off 
time Ih for scale-free networks with BA parameter m = 3 under different 
external fields, as indicated. 

external fields and observables: (a) the total magnetization, (b) the absolute 
magnetization, and (c) the fraction of final states with a majority of up spins. 
The results shown correspond to scale-free networks with BA parameter m = 
3. From a qualitative point of view, these results compare well to the small- 
world data shown above (see Figure 6 and discussion). 

Finally, let us consider in Figure 11 the ratio fupiT = 1)/ fupiT = 4) as a 
function of tn for different fields, focusing again on scale-free networks with 
BA parameter m = 3. For H = 0.1, the ratio shows a slight increase for small 
values of up to about — S x 10~^, while it decreases monotonically for 
higher values of tH and reaches a minimum at t// = 1 (horizontal arrow). This 
behavior resembles the observations made above for SWNs (recall Figure 
7), although the effect apparent here is milder. For H = 0.5, the ratio 
peaks at tn ~ 1-5 x 10^^ (vertical arrow) and then decreases for larger 
values oHh- This confirms that, irrespective of the topological details of the 
complex network substrate, short, intense campaigns are the best strategy to 
successfully convey a message across a highly cohesive society. Incidentally, 
notice also that, unlike the SWN case, the 1 tail seems to be roughly 

independent of the external field. 

at Tc « 2 [35], the chosen temperatures roughly correspond to T = O.STc and T = 2Tc, 
respectively. 
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4 Conclusions and Outlook 



In order to investigate advertising and irreversible opinion spreading phenom- 
ena, we studied the magnetic Eden model (MEM) growing on small-world and 
scale-free networks with externally applied magnetic fields. In this nonequi- 
librium model, the opinion or decision of an individual is affected by those 
of their acquaintances, but opinion changes (analogous to spin flips in an 
Ising-like model) do not occur. 

The model studied in this work could be realistically applied to socio- 
logical scenarios in which individuals are subject to highly polarized, short 
term, binary choice situations under the influence of advertising/marketing 
mass media campaigns. Given these conditions, opinions are not expected 
to fluctuate and "thermalize" . Some examples are binary voting situations, 
such as a ballotage or referendum, and market competition scenarios between 
two leading products, where the validity of this model could be empirically 
tested using polls, surveys and/or sales performance data. 

The interplay and competition between the intrinsic local disorder and 
the externally applied ordering field was investigated by means of order- 
parameter probability distributions and ensemble averages, leading to a tem- 
perature vs. magnetic field order- disorder phase diagram. The small-world 
network's shortcut fraction was found to shift the effective temperature of 
the system, suggesting a simple scahng relation that allowed a collapse of all 
small-world network data. Later, the effects of advertising campaigns with 
variable duration was discussed. The best cost-effective strategy to achieve 
a desired majority consensus was found to depend crucially on the intrin- 
sic cohesion tendency of the society. In the scenario of strongly cohesive 
social groups, intense, short campaigns are observed to be more successful 
than mild, long ones. However, the opposite conclusion holds in the case of 
weakly cohesive societies. 

The observed phenomena, first presented and discussed for small-world 
network substrates, were compared to corresponding results obtained using 
scale-free network geometries. We found that, irrespective of the details of the 
underlying complex network topology, consensus and order can be reached for 
any level of intrinsic cohesion and interconnectedness. Moreover, the effects 
of advertising campaigns with variable duration, as well as our conclusions 
regarding optimal strategies to achieve a majority consensus, remain valid 
irrespective of the specific type of complex network substrate. 

The magnetic Eden model, being different in nature from related equilib- 
rium spin systems such as the Ising model, can certainly provide valuable, 
complementary insight into dynamical and critical aspects of the spreading 
of opinions in a society. Since a characteristic feature of the magnetic Eden 
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model is the absence of thermalization, it may be interesting to study a 
hybrid model with an external parameter controlling the rate of thermaliza- 
tion relative to the rate of growth. This generalized model could be tuned 
from pure MEM-like growth to pure Ising-like thermalization, which might 
be an appropriate representation of sociological scenarios where the opinion 
of individuals is neither completely frozen nor fully thermalized. 

Hopefully, the present findings will thus contribute to the growing in- 
terdisciplinary efforts in the fields of sociophysics, complex networks, and 
nonequilibrium statistical physics, and stimulate further work. 
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